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Background 

Acoustic  rays  or  modes  propagating  in  the  shallow  water  waveguide  combine  constructively  and 
destructively  to  produce  richly  structured  interference  patterns.  In  the  early  1980s,  Russian 
scientists  showed  how  the  primary  features  of  these  interference  patterns  could  be  distilled  into  a 
single  scalar  parameter,  the  so-called  waveguide  invariant.  In  recent  year,  the  waveguide 
invariant  concept  has  become  central  to  many  signal  processing  schemes  for  both  passive  and 
active  sonar.  A  difficulty  in  applying  the  concept  is  that  one  must  know  the  waveguide 
invariant’s  numerical  value.  The  canonical  value  in  shallow  water  is  1.0,  but  several  experiments 
have  shown  that  a  larger  number  is  often  appropriate  when  the  sound  speed  profile  has  a  sharp 
gradient.  Deviations  from  the  canonical  value  become  more  pronounced  as  the  frequency  is 
increased.  These  deviations  are  of  more  than  academic  interest:  errors  in  the  assumed  numerical 
value  for  the  waveguide  invariant  translate  directly  into  errors  in  the  estimated  range  to  a  target. 

Work  Completed 

The  waveguide  invariant  has  traditionally  been  regarded  as  a  low  frequency  technique.  Under 
current  support,  the  objective  has  been  to  address  issues  relevant  to  extending  the  method  into 
the  mid- frequency  (1-10  kHz)  regime.  The  first  objective  was  to  develop  an  algorithm  for 
estimating  the  numerical  value  of  the  waveguide  invariant  directly  from  acoustic  data  without 
having  a  priori  knowledge  of  a  source’s  range  or  depth.  The  second  objective  was  to  study  how 
signal  fluctuations  (scintillation)  caused  by  random  variability  in  the  water  column  from  such 
features  as  internal  waves  would  affect  the  waveguide  invariant  concept.  The  work  has  made  use 
of  existing  data  sets  collected  under  ONR  support:  the  Shallow  Water  2006  (SW06)  experiment 
and  the  2011  Gulf  Experiment  (GulfExl  1) 

As  listed  below,  the  completed  work  has  been  documented  in  7  journal  publications  and  5  invited 
conference  presentations. 

Research  Highlight 

Existing  algorithms  allow  one  to  estimate  only  the  ratio  between  the  invariant  and  range:  to  get 
the  source’s  range,  one  must  know  the  numerical  value  of  the  invariant.  The  primary  result 
produced  under  current  funding  is  a  new  algorithm  for  isolating  range  from  the  waveguide 
invariant  thereby  permitting  one  to  calculate  both  quantities  and  not  just  the  ratio.  The  algorithm 
works  by  twice  beamforming  the  measured  acoustic  field.  Beamforming  conventionally,  the 
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processor  first  focuses  at  the  true  bearing  of  the  source.  The  field  is  then  beamformed  with  each 
hydrophone  in  the  array  evaluated  at  a  slightly  different  frequency.  The  frequency-shifted  output 
will  shift  the  processor’s  focus  to  a  different  bearing.  The  extent  to  which  the  focus  shifts 
between  the  two  calculations  provides  a  direct  measurement  of  the  waveguide  invariant  that  is 
independent  of  the  source’s  range.  With  this  independent  knowledge  of  the  waveguide  invariant, 
one  can  then  calculate  the  range  unambiguously. 

The  method  of  twice-beamforming  the  field  was  first  demonstrated  successfully  with  numerical 
simulations.  It  was  further  shown  that  the  method  could  provide  a  coarse  estimate  of  source’s 
depth  as  it  could  distinguish  near-surface  sources  from  sources  that  were  located  below  the 
mixed  layer.  Analysis  of  data  collected  during  GulfExll  provided  additional  support  for  the 
proposed  algorithm. 

The  complete  paper  is  included  as  part  of  this  report. 

Professional  Activities  and  Honors 

Acoustical  Society  of  America.  Technical  Committee  service:  Underwater  Acoustics  (three 
tenns,  1998-2004  and  2009-2012),  Acoustical  Oceanography  (four  terms,  2002-2014).  Technical 
Program  Organizing  Meeting  (TPOM)  representative  for  Underwater  Acoustics  and  Acoustical 
Oceanography,  Hong  Kon.  Organizer  and  Co-chair  of  Special  Session  “Underwater  acoustic 
communications  and  networking,”  Joint  Acoustical  Society  of  America/Chinese  Acoustical 
Society  Meeting,  Hong  Kong  May  14-18,  2012. 

Excellent  Reviewer  Award,  IEEE  Journal  of  Oceanic  Engineering,  2012. 

Invited  lecture  series,  Xiamen  University,  Xiamen,  China,  November  27-30,  2011. 

Invited  lecture  series,  Ocean  University,  Qingdao,  China,  November  22-26,  2011. 

Invited  lecture,  Harbin  Engineering  University,  Harbin,  China,  July  19,  2010. 

Invited  lecture  series,  Zhejiang  University,  Hangzhou,  China,  July  12-14,  2010. 
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Abstract:  The  waveguide  invariant  summarizes  the  pattern  of  con¬ 
structive  and  destructive  interference  between  acoustic  modes  propagat¬ 
ing  in  the  ocean  waveguide.  For  many  sonar  signal-processing  schemes, 
it  is  essential  to  know  the  correct  numerical  value  for  the  waveguide 
invariant.  While  conventional  beamforming  can  estimate  the  ratio 
between  the  waveguide  invariant  and  the  range  to  the  source,  it  cannot 
unambiguously  separate  the  two  tenns.  In  the  present  work,  striation- 
based  beamforming  is  developed.  It  is  shown  that  the  striation-based 
beamformer  can  be  used  to  produce  an  estimate  for  the  waveguide  invar¬ 
iant  that  is  independent  of  the  range.  Simulation  results  are  presented. 
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1.  Introduction 

Consider  an  underwater  passive  sonar  scenario  where  a  horizontal  array  measures  the 
acoustic  field  produced  by  a  distant  source.  Mapped  as  a  function  of  frequency  m  and 
range  r,  the  observed  acoustic  intensity  /  often  exhibits  striations.  These  striations, 
alternating  bands  of  high  and  low  intensity,  result  from  constructive  and  destructive  in¬ 
terference  between  the  propagating  acoustic  modes  supported  by  the  ocean  waveguide. 
Chuprov1  developed  a  simple  formula  for  the  slope  dcoldr  of  the  striations: 


dm 

dr 


—di/dr 

di/dm 


mfi/r, 


(1) 


where  jl  is  the  so-called  waveguide  invariant.  Typically,  one  would  know  the  frequency 
and  measure  the  slope  allowing  the  fl/r  ratio  to  be  calculated. 

While  the  ratio  between  the  waveguide  invariant  and  the  range  may  be  readily 
determined,  isolating  the  range  from  the  numerical  value  for  />  is  more  problematic.  In 
practice,  one  might  simply  assume  that  /I  is  known  and  then  calculate  r  from  the  jilr 
ratio.  The  canonical  value  in  shallow  water  is  j8=  1,  but  this  choice  may  prove  far  too 
coarse.  If  both  the  source  and  array  are  below  the  thermocline,  both  simulations2  and 
experiments3’4  have  shown  that  /?=  1.3  or  1.4  is  more  appropriate.  Erroneously  assum¬ 
ing  /I  =  1  when  it  is  actually  /?=  1.3  translates  directly  into  a  30%  error  in  the  estimated 
value  for  the  range. 

In  the  present  work,  a  method  is  proposed  for  isolating  the  waveguide  invari¬ 
ant  and  producing  an  estimate  for  its  numerical  value  that  is  independent  of  the  range. 
The  method  relies  on  twice  beamforming  the  acoustic  pressure  measured  on  a  horizon¬ 
tal  array  due  to  a  distant  source.  The  pressure  field  is  first  beamformed  conventionally 
at  a  single  frequency  to  achieve  a  focus  at  the  correct  bearing  of  the  source.  The  field 
is  then  beamformed  again  but  with  each  point  on  the  array  evaluated  at  a  slightly 


EL76  J.  Acoust.  Soc.  Am.  130  (2),  August  201 1 


©  201 1  Acoustical  Society  of  America 


Distribution  Statement  A:  Approved  for  public  release;  distribution  is  unlimited. 


Downloaded  25  Jul  2011  to  128.95.76.220.  Redistribution  subject  to  ASA  license  or  copyright;  see  http://asadl.org/journals/doc/ASALIB-home/info/terms.jsp 


D.  Rouseff  and  L  M.  Zurk:  JASA  Express  Letters 


[DOI:  10.1121/1.3606571] 


Published  Online  22  July  201 1 


different  frequency.  The  frequency  offset  is  selected  based  on  the  observed  slope  of  the 
intensity  striations.  It  is  shown  that  the  focus  will  shift  from  what  is  observed  with  con¬ 
ventional  beamforming  and  that  the  extent  of  the  shift  gives  a  measure  of  />  that  is  in¬ 
dependent  of  r.  The  method  also  gives  an  indirect  coarse  estimate  of  the  source’s 
depth. 

To  begin  the  development,  assume  a  point  source  at  r0  =  (r0  cos  </)0, 
r0  sin  cj)0,z0)  in  the  ocean  waveguide.  The  resulting  field  is  measured  on  a  horizontal 
line  array  of  length  L  oriented  along  the  >-axis,  —  L/2  <  y  <  L/2,  at  depth  za.  Assum¬ 
ing  a  range-independent  environment  and  an  exp(— imt)  convention,  the  measured 
complex  field  is  expanded  in  normal  modes  as 

p(y,  ffl)  =  exp(i^mr),  (2) 


where  each  mode  'P„,  has  corresponding  horizontal  wavenumber  and  certain  unim¬ 
portant  scaling  factors  are  neglected.  It  is  assumed  that  the  source  is  coherent  over  the 
narrow  frequency  band  of  interest.  The  range  r  to  a  point  on  the  array  is 

r  =  [(ro  cos  4 >0)2  +  [y  -  r0  sin  0O)2]1/2 
~r0-y  sin  4>0  +  y2  /  (2r0)  (3) 


to  the  term  second  order  in  y  that  expresses  wavefront  curvature. 

Since  /  =  |p|2,  the  slope  of  the  observed  intensity  follows  from  Eqs.  (1 )-  (3): 

dm  dm  dr 

7F=*S^“(/!/r")sm^  <4) 

where  the  distant  source  is  assumed  not  too  near  broadside,  4>o  =  if  so  the  wavefront 
curvature  term  can  be  neglected.  Solving  the  separable  differential  equation  yields 

In  (co/too)  =  ~(P/ro)  sin  {4>0)y,  (5) 

that  may  be  approximated  for  frequencies  near  the  center  frequency  m0  as 

m  =  m0  —  Sy,  where  S  =  m0(P/r0)  sin  cj)0.  (6) 


Equation  (6)  says  that  the  level  curves  of  constant  intensity  as  observed  along  the  array 
will  map  onto  straight  lines  of  slope  S.  Consequently,  incoherent  processing  using  only 
the  intensity  is  sufficient  to  estimate  S.  The  numerical  value  for  S  might  be  calculated 
by  spectral  methods,  as  originally  suggested  by  Chuprov,1  or  perhaps  by  using  Radon 
transforms.5  Regardless  of  the  specific  method  used  to  do  the  calculation,  the  remain¬ 
der  of  this  report  assumes  that  the  value  for  S  is  known. 

Having  first  processed  the  measured  field  incoherently  to  deduce  S,  it  is  now 
processed  coherently  by  beamforming.  Define  the  beamformer  output  in  look  direction 
<p  as 


B{4>)  =  \b(4>)\>  with  b{4>)  =  L  1 


■L/2 

pLv,co)w*(y,co;4>)dy. 

-L/2 


(7) 


For  a  conventional  Bartlett  beamformer,  the  weights  are 

w*{y,  m;  4>)  =  exp[z'(m/c0)  sin (<p)y\,  (8) 


where  c0  is  the  reference  sound  speed.  Substituting  Eqs.  (2)  and  (8)  into  (7)  and 
rearranging  terms  yields 
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b{4>)  =  X  CmDm{w,4>)- 

m 


(9) 


The  terms  slowly  varying  in  y  and  to  are  lumped  into  Cm.  The  rapidly  varying  phase 
terms  are  included  in  Dm\ 


Dm{co;  <p)  =  L  1 


■L/2 

exp \i^m((o)r+i{(o/co)  sin(<j!>)y]  dy, 

—L/2 


(10) 


where  the  frequency  dependence  of  the  horizontal  wavenumber  cjm  has  been  made 
explicit. 

In  conventional  beamforming,  the  frequency  is  fixed  at  some  constant  value, 
say,  co  =  co0.  In  the  present  work,  a  striation-based  beamformer  is  considered  where 
different  points  along  the  array  are  evaluated  at  different  frequencies.  The  frequency 
varies  linearly  across  the  array  based  on  the  observed  slope  S  of  the  intensity  stria- 
tions.  To  maintain  generality,  let  the  frequency  across  the  array  be 

co  =  ct>o  —  oiSy,  (11) 


where  a  is  an  adjustable  parameter.  For  conventional  beamforming,  a  =  0  and  the  fre¬ 
quency  is  constant  across  the  array.  For  beamforming  along  the  striations,  a=l;  see 
Eq.  (6).  The  frequency-dependent  modal  wavenumber  across  the  array  can  be  approxi¬ 
mated  by  a  Taylor’s  series  expansion.  To  first  order, 

o  -  aSy)  ~  <L(® o)  +  (d£m/d(o)(-<xSy) 

=  co0/vm-  (oiS/um)y,  (12) 


where  vm  and  um  are  the  phase  and  group  velocities,  respectively,  for  mode  m  at 
co  =  coq.  The  expansion  in  Eq.  (12)  should  retain  the  quadratic  term  to  be  strictly  con¬ 
sistent  with  the  expansion  in  Eq.  (3).  The  quadratic  term  in  expansions  like  Eq.  (12)  is 
usually  neglected;  see  the  discussion  in  Ref.  6  for  the  limitation  associated  with  this 
approximation. 

Combining  the  above  equations  yields  for  the  striation-based  beamformer 


Dm(co0  -  ctSy;  </>)  =  L  1 


•L/2 

exp  [/(d0  +  diy  +  d2y2)]  dy, 

-L/2 


(13) 


where  terms  have  been  grouped  in  like  order  of  y  with  coefficients 

do  =  (fflo/v„,)7o, 


(14) 


di  =  (coo/co)  sin (0)  -  (®o/vm)  sin(<j f>0)  -  a Sr0/um 

=  (mo/co)  sin(0)  -  {co0/vm)  sin(0o)[l  +  a P(vm/um)} ,  (15) 

d2  =  (®o/vm)/(2r0)  +  (aS/M,„)[sin((/)0)  -  (um/c0)  sin (</>)].  (16) 


It  is  instructive  to  consider  the  linear  coefficient  d|  in  greater  detail.  In  transitioning 
between  the  two  forms  of  Eq.  (15),  the  definition  of  the  striation  slope  S,  Eq.  (6),  has 
been  applied.  Since  S  is  proportional  to  the  ratio  /i/r0,  the  product  Sr0  appearing  in 
Eq.  (15)  is  independent  of  r0.  The  key  point  is  that  d|  therefore  depends  on  the  wave¬ 
guide  invariant  but  is  independent  of  range;  /i  has  been  isolated  from  r0  by  striation- 
based  beamforming. 

Equation  (13)  can  be  evaluated  exactly  in  terms  of  Fresnel  integrals.7  To  gain 
greater  insight  into  the  result,  it  is  useful  to  assume  that  the  source  is  in  the  far  field  of 
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the  receiving  array.  In  the  far  field,  the  quadratic  term  in  Eq.  (13)  is  neglected  and  the 
integral  reduces  to 

Dm(a>o  —  uSy\  (p)  =  e'e"sinc(0iL/2),  (17) 


where  sinc(x)  =  sin(jc) /x.  The  magnitude  is  a  maximum  when  d\  =  0  implying  focusing 
when  the  beamformer  look  direction  (p  satisfies 

sin (<p)  =  (co/v,„)  sin(0o)[l  +  aj8(vm/wm)] 

«  sin(0o)(l  +  «/?),  (18) 

where  the  approximate  form  ignores  the  small  differences  between  c0  and  the  phase 
and  group  velocities  to  yield  a  result  that  is  independent  of  the  mode  index. 

Equations  (6)  and  (18)  suggest  a  three  step  processing  strategy  for  estimating 
the  waveguide  invariant.  First,  process  the  field  incoherently  to  estimate  the  slope  5  of 
the  intensity  striations.  Second,  beamform  conventionally  with  a  =  0.  The  beamformer 
will  focus  at  the  true  bearing  of  the  source,  cp  =  cpQ.  Finally,  perform  striation-based 
beamforming  with  a  /  0.  The  focus  will  shift  to  some  new  (p  /  (pf)  but,  since  the  true 
<p0  has  already  been  estimated,  the  location  of  the  new  focus  can  be  used  in  Eq.  (18)  to 
solve  for  /?.  The  value  for  the  waveguide  invariant  is  thereby  estimated  without  know¬ 
ing  the  range  or  depth  of  the  source. 

The  numerical  value  for  the  adjustable  parameter  a  should  be  selected  with 
some  care.  The  choice  a  =  1  is  attractive  because  it  means  a  strong  signal  across  the 
array  without  fades.  This  should  improve  performance  in  a  noisy  environment.  When 
the  source  is  at  a  bearing  too  far  from  array  broadside,  however,  the  choice  a  =  1  will 
fail  to  produce  a  focus.  For  example,  no  beamformer  look  direction  cp  will  satisfy  Eq. 
(18)  if  a.  =  ft  =  1.0  and  \cp0\  >30°.  If  conventional  beamforming  (a  =  0)  reveals  a  source 
far  removed  from  broadside,  the  striation-based  beamformer  should  use  a  value  for  a 
that  ensures  that  Eq.  (18)  will  have  a  solution.  In  addition,  a  must  be  selected  so  that 
the  total  frequency  shift  across  the  array,  ocSL,  is  within  the  coherence  bandwidth  of 
the  source. 

A  specific  numerical  simulation  serves  to  illustrate  the  striation-based  beam¬ 
forming  concept.  Consider  an  environment  typical  of  the  East  China  Sea  as  observed 
during  a  2001  experiment.8  The  105  m  deep  water  column  has  a  40  m  deep  surface 
mixed  layer  above  a  sharp  thermocline.  The  sound  speed  contrast  between  the  surface 
mixed  layer  and  bottom  of  the  water  column  is  12  m/s.  Assume  a  128  element  horizon¬ 
tal  receiving  array  with  0.5  m  spacing  (L  =  63.5  m)  positioned  at  depth  za  =  60  m.  The 
source  is  positioned  at  range  r0=  5  km  and  bearing  relative  to  array  broadside  cp0 
=  20°.  The  initial  source  depth  is  zo  =  60  m  so  that  both  the  source  and  the  receiving 
array  are  positioned  below  the  part  of  the  thermocline  with  a  sharp  gradient.  With 
both  the  source  and  the  receiving  array  below  the  main  thermocline,  the  low-order 
acoustic  modes  that  have  an  upper  turning  depth  below  the  sea  surface  will  contribute 
strongly  to  the  acoustic  field.  In  such  a  scenario,  previous  results2-4  show  that  /?  >  1.0 
might  be  expected. 

Figure  1  shows  the  simulated  intensity  as  calculated  along  the  array  over  a  50 
Hz  band.  The  calculation  was  made  using  a  normal  mode  code  and  the  plot  has  20  dB 
dynamic  range.  Striations  are  plainly  evident.  A  plausible  estimate  for  the  slope  of  the 
striations  is  5=0.089  Hz/m  implying  a  SL  =  5.7  Hz  frequency  shift  over  the  length  of 
the  array.  Superimposed  on  the  image  are  lines  corresponding  to  a  =  0  and  a  =  1  using 
reference  frequency  mJ2n=  1  kHz  in  Eq.  (11). 

A  discrete  Bartlett  beamformer  was  implemented  and  Fig.  2  shows  the  results. 
Beamforming  conventionally  (oe  =  0),  the  figure  shows  a  focus  at  look  direction 
<p  =  (po  =  20°  as  expected.  The  plot  is  normalized  so  the  conventional  B(cp0)=  1.0.  The 
field  was  then  beamformed  along  the  striation  (a  =  1 )  with  the  same  normalization 
applied  and  the  result  plotted.  As  expected,  the  output  at  the  focus  is  higher  for  the 
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Fig.  1 .  (Color  online)  Acoustic  intensity  observed  across  horizontal  array  as  function  of  frequency.  Superim¬ 
posed  lines  a  =  0  and  a  =  1  denote  trajectories  for  conventional  and  striation-based  beamforming,  respectively. 
Dynamic  range  is  20  dB.  See  text  for  simulation  details. 


-30  -20  -10  0  10  20  30 

distance  along  array  (m) 


striation-based  beamformer.  More  importantly,  the  location  of  the  focus  has  shifted  to 
<£  =  52°.  Using  (j>  =  52°  and  </>0  =  20°  in  Eq.  (18)  yields  a  value  for  the  waveguide  invar¬ 
iant  £  =  1.3,  a  value  consistent  with  experimental  results3'4  when  both  source  and  re¬ 
ceiver  are  below  the  main  thermocline. 

The  simulation  was  repeated  with  all  the  same  parameters  except  the  source 
depth  was  changed  to  z0  =  10  m.  With  the  source  in  the  surface  mixed  layer,  the  low- 
order  modes  that  cause  /?  >1.0  are  not  excited.  Striations  are  again  evident,  but  with  a 
shallower  slope  than  is  observed  in  Fig.  1;  S=  0.069  Hz/m  is  a  reasonable  value.  The 
conventional  and  striation-based  beamformers  focus  at  <£  =  <£o  =  20°  and  <£  =  43°, 
respectively.  Using  these  angles  in  Eq.  (18)  yields  £=1.0  for  the  shallow  source. 

Because  the  striation-based  beamformer  is  able  to  quantify  £,  it  also  provides 
a  coarse  estimate  of  the  source  depth,  at  least  to  the  extent  that  it  can  distinguish 
between  shallow  and  deep  sources.  A  caveat  is  that  the  present  analysis  has  assumed  a 
range-independent  environment.  Range  dependence  from  such  environmental  factors 
as  internal  waves  will  cause  the  acoustic  modes  to  couple7  and  make  source  depth  esti¬ 
mation  more  problematic. 

To  summarize,  while  conventional  beamforming  can  determine  a  source’s 
bearing  and  the  £/r0  ratio,  it  is  unable  to  separate  the  waveguide  invariant  £  from  the 
range  r0.  Striation-based  beamforming  breaks  the  £/r0  ambiguity  and  provides  an  inde¬ 
pendent  estimate  for  the  waveguide  invariant. 
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Fig.  2.  (Color  online)  Beamforming  results.  True  source  at  bearing  0O  =  20°.  Shift  in  location  of  focus  from 
conventional  (a  =  0)  to  striation-based  (a  =  1 )  beamforming  can  be  used  to  estimate  the  waveguide  invariant 
using  Eq.  (18).  See  text  for  details. 
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